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Abstract — Linear receivers offer a low complexity option for 
multi-antenna communication systems. Therefore, understanding 
the outage behavior of the corresponding SINR is important in 
a fading mobile environment. In this paper we introduce a large 
deviations method, valid nominally for a large number M of 
antennas, which provides the probability density of the SINR of 
Gaussian channel MIMO Minimum Mean Square Error (MMSE) 
and zero-forcing (ZF) receivers, with arbitrary transmission 
power profiles and in the presence of receiver antenna correla- 
tions. This approach extends the Gaussian approximation of the 
SINR, valid for large M asymptotically close to the center of the 
distribution, obtaining the non-Gaussian tails of the distribution. 
Our methodology allows us to calculate the SINR distribution 
to next-to-leading order (0(1/M)) and showcase the deviations 
from approximations that have appeared in the literature (e.g. 
the Gaussian or the generalized Gamma distribution). We also 
analytically evaluate the outage probability, as well as the 
uncoded bit-error-rate. We find that our approximation is quite 
accurate even for the smallest antenna arrays (2 x 2). 

Index Terms — Gaussian approximation, information capacity, 
large-system limit, multiple-input multiple-output (MIMO) chan- 
nels. 

I. Introduction 

Multi-antenna systems have been known (TJ, J2] to offer 
considerable advantages, not only at the link-level, providing 
higher multiplexing gains and increased robustness through di- 
versity, but also at a system-level by allowing a more effective 
interference mitigation in a multi-user setting. It is therefore no 
surprise that next generation wireless communications systems 
will include multi-antenna systems O in order to capitalize 
on these benefits. To obtain the full advantages from multiple 
antennas, it is necessary to have an optimal receiver structure, 
which however is quite complex to implement in real systems. 
Instead, low complexity, albeit suboptimal, linear receivers 
offer a practical alternative. 

Such receivers include the so-called MMSE (minimum 
mean square error) and the zero-forcing (ZF) receivers, as 
well as a new class of receivers recently proposed [4] called 
moment-based receivers. In addition to the simplification due 
to the linearization of the received signal operation, the re- 
ceived signal may then be iteratively treated to cancel the 
interference from other antennas. However, in many cases, 
even this may impose significant complexity. An even simpler 
receiver structure can be constructed, in which, after the linear 
spatial equalization the data is decoded in a single-input single- 
output fashion Q, j6). Here we will focus on the latter, 
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especially since we are interested in the cellular context, with 
separated transmitter antenna arrays in the uplink with a multi- 
antenna receiver terminal. 

The throughput performance depends on the ability of the 
linear receiver structure to mitigate interference. One very 
useful method to quantify the performance is through the 
asymptotic analysis of the signal to interference and noise ratio 
(SINR) for the receiver in the limit of large antenna numbers 
using tools from random matrix theory. Its application was 
initially spearheaded in the context of Direct-Sequence Code- 
Division-Multiple-Access (DS-CDMA) where the effective 
channel consists of the matrix of pseudorandom codes. In this 
direction, the first breakthrough was made by Q, 0, who 
showed that in the infinite matrix-size limit the SINR of a 
fixed random channel realization converges to its mean. Later, 
similar results were obtained for more general channels (9), 
[ 1 1 . More recently, the effectiveness of linear receivers were 
analyzed in terms of the total throughput from all transmitting 
nodes in the asymptotic limit lfT71 - lfT3l . 

Nevertheless, one often needs to assume that the fading 
channel is "quasi-static", i.e. varies in time much more slowly 
than the typical coding delay. In this case the channel matrix 
and hence the SINR have to be considered as random quanti- 
ties. In this regime, the relevant performance metric is the "rate 
(or SINR) versus outage probability" tradeoff lfl4ll . captured by 
the cumulative distribution function of the SINR. This situation 
is especially relevant in the context of multi-antenna channels, 
when the number of antennas is usually much smaller than the 
size of the CDMA codes. 

In a seminal work |L5l the authors proved the asymptotic 
normality of the SINR for the MMSE and ZF receivers when 
all transmitters have equal power. The normality of the SINR 
was later extended to the normality of the multiple access inter- 
ference (MAI) of CDMA channels lfl6ll and a variety of linear 
receivers 1171 . More recently, ifTsl , 11191 showed the normality 
of the MMSE SINR, including the case of the mismatched 
receiver. Interestingly, [18] showed also that the logarithm 
of the SINR becomes asymptotically normal. Unfortunately, 
and in contrast to the total mutual information, the Gaussian 
approximation for the SINR is extremely inaccurate, unless 
the number of antennas is quite large. As a result, inspired 
by the fact that the SINR for the equal power MIMO ZF 
receiver has a Gamma distribution 1151 . (20 1, several works 
were devoted to approximating the SINR statistics with other 
distributions, such as the Beta distribution for the SINR of the 
CDMA ZF receiver lETI . 11221 . or the Gamma and generalized 
Gamma distributions |23|-[26|, in which case their first three 
moments were fitted to match the actual SINR distribution. 
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Nevertheless, this methodology, although perhaps providing 
good agreement under certain conditions, is ad-hoc and does 
not offer any intuition on the SINR statistics. The same can 
be argued for the calculation of the exact probability density 
function (PDF) and the cumulative density function (CDF) of 
the MMSE SINR ll27l using ratios of determinants, a method 
however which is only valid for uncorrelated channels at the 
receiver. 

In this paper, we take a different approach. Instead of trying 
to prove Gaussian behavior close to the peak of the distribution 
of SINR, we develop a large-deviations methodology, which 
allows us to calculate the distribution of the SINR arbitrarily 
far from its most probable, mean value. The success of 
our method lies on the fact that we can exactly express 
the moment generating function (MGF) of the SINR as the 
moment generating function of the difference of two correlated 
MIMO mutual information functions. Taking advantage of 
the robustness of the Gaussian approximation of the MIMO 
mutual information we obtained an expression of the MGF 
of the SINR correct to 0(1/ M). We are then able to obtain 
the full distribution of the SINR for both MMSE and ZF with 
similar precision. It is therefore no surprise that our results 
are very close to the exact ones down to the smallest MIMO 
systems (2 x 2). It is worth mentioning a related recent work 
1 28 1 in which we used the Coulomb Gas method ||29l to 
calculate the leading term 0(M) in the exponent of the SINR 
for uncorrelated channels. In that work we demonstrate that 
the large deviations tails are determined by the behavior of a 
single singular value of the channel matrix. 

Outline: In the next section we present the channel model 
and introduce the MMSE and the ZF SINR. In Section [111] 
we present our analytical results, providing the PDF, CDF 
and BER for both MMSE and ZF SINRs. In Section [TV] we 
demonstrate their validity numerically and we conclude in 
Section [V] Appendices [A] [B] [C] and [0 contain details on the 
proofs of Lemma [TJ Propositions [JJ [2] and [3] respectively. 



is an M -dimensional vector with i.i.d. entries ~ CAf(0, 1). 
Ro is the A/-dimensional receive-side correlation matrix of 
the channel originating from user 0, normalized so that 
trRo = M. Similarly, the channel matrix from the kth 

1 /2 

user is H k Gfc, where Gfe is a M x Nk matrix with i.i.d. 
elements ~ CJ\f(0, 1/Nk) and R& has the same interpretation 
and properties as Ro. To be concrete, we will assume that 
all correlation matrices R&, for fc = 1,...,K are positive 
semidefinite, while Ro is positive definite. Also, we assume 
that their eigenvalue spectra converge to proper probability 
distributions for large A/. We will be interested in calculating 
the SINR of transmitter in the presence of the other trans- 
mitters and noise. For notational convenience we also define 
the matrix Ho = [Rj^ 2 gi y/pT, . . . , R^ 2 Gif \[pk\ as well as 
the matrix H = [R^go^/po", H ]. 

This channel model describes a set of transmitting antennas 
dispersed in a cellular setting with their signal arriving pos- 
sibly from different mean angles and/or with different angle- 
spreads at the receiver array, thereby having different receive 
correlation matrices. Of course, not all correlation matrices 
need to be different, e.g. if some of the interfering antennas 
are collocated. To obtain analytic results we will take the limit 
of large M and Nk (k = 1 , . . . , K), with the ratios 



as well as the number of arrays K fixed in that limit. In the re- 
mainder of the paper the term "large M limit" will denote both 
Nk and M going to infinity, while keeping the corresponding 
ratios rik constant and finite. For notational convenience, we 
define iV tof = X^tLo ^h- Despite the assumptions above, we 
will apply and test our results in the case when M and K are 
not too large and Nk = 1. 



II. Channel Model 

In this section we define the channel model. The receiver 
array has M antennas, receiving the signal from K + 1 
transmitter arrays, not necessarily collocated. Without loss of 
generalitjQ we assume that the 0th transmitter has a single 
antenna, while the remaining K transmitters have Nk antennas 
each for k = 1, . . . , K. The M-dimensional received signal 
vector y can be written as 

K 

y = Ro /2 g x + Rfe /2 G fc x fc + z (1) 

fe=i 

In the above equation z is the noise vector, with complex 
Gaussian elements ~ CJ\f(0, 1). The transmitted signal am- 
plitudes xq and Xfc have i.i.d. elements with variance po 
and PfcljVfc respectively, where pk are the average transmitted 
power per antenna from the fcth array with k = 1,...,K. 

1 /2 

The channel vector from transmitter is R go, where go 

1 For example, we may assume that the zeroth and the first arrays correspond 
to the same transmitter 



A. MMSE Receiver 

The SINR of the 0-th MMSE transmitter above can be 
expressed as 

7(H) = ^So R o L R o So (3) 

L = I M + H H (4) 

with the second line serving as the definition of L. Our 
objective is to evaluate the probability density function of 
7(Hq), omitting the Ho dependence when obvious. 



B. ZF Receiver 

The SINR of the zero-forcing (ZF) receiver can be obtained 
in a similar fashion. In this case, we focus only in the case 
M > N to t- Then the SINR for this receiver can be expressed 
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as a limit of the standard MMSE SINR (0 as follows 

Po 1 



7(H) = 



M 
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The inverse of the matrix in the right-hand side of the first 
equality is finite only for M > N to t with probability one. The 
second equality results from taking the z —> + limit. The 
third equality above results from the matrix inversion lemma 
[ 30 1 . Following the same argumentation as in Section Hi- Al we 
obtain the moment generating function as in ©, with H Hj in 
(0 replaced by z HqH-o- The expression in the fourth line, 
easily derived from the third, showcases the singular nature of 
the z — > limit, which focuses on the projection of the kernel 
of HoHj to go, which, for M > N tot is guaranteed to be 
non-empty. For compactness, below we will continue to use 
this dummy variable z, setting it equal to z — 1 and z = + 
for the cases of the MMSE and ZF SINR, respectively. 

III. Results 

In this section we will go through the basic steps of the 
calculation of the probability distribution (PDF), the outage 
distribution (CDF) and the BER of the SINR denoted by 7. We 
start with a very useful first result for the moment generating 
function of 7. 

Lemma 1 (MGF of 7). The moment generating function of ^ 
for the MMSE (0 and the ZF case (0 can be written in the 
following form 

9m(s) = E H o{e- A/(s ' Ho) } (6) 
where A/(s, Hq) is given by 



A/(s,H ) = lndet 



In det 



z(I M + spoRo) + H H 



zl 



M 



H H 

The parameter z takes the z = 1 for the MMSE SINR (0 and 
the limiting value z = + for the ZF SINR (0. 

Proof: See Appendix lAl 
Remark 1. Once again we see that the limit z — >• + in (0 is 
not trivial because the matrix HoHq has a non-empty kernel. 

Remark 2. The usefulness of this result is that it makes the 
connection of the moment generating function of the SINR to 
a difference of mutual information functions for the remaining 
K users. This will allow us to take advantage of the Gaussian 
behavior of this difference of mutual informations [31 1, 11321 
close to their ergodic values, in order to analyze the large 
deviations of the distribution of 7 arbitrarily far away from its 
ergodic value. Note that the above argument holds for general 
Ho, as long as the logdets difference above remains Gaussian, 
as e.g. in 



(7) 



A. Derivation of PDF 

We will now obtain the probability distribution density of 
the SINR. This density may be expressed as an expectation of 
a Dirac <5-function as follows: 

-1 



P M {l) = E 



H 



7' 



z P0i 



zl 



M 



H H 



go 



(8) 



The parameter z will take the value of z = 1 for the case 
of the MMSE SINR introduced in Section III-A1 while, as 
discussed in Section ITl-BI the z = + limit will correspond to 
the ZF SINR. The following proposition provides an analytic 
expression of the probability density of the SINR, valid for all 
7 > in the large M limit. 



Proposition 1 (PDF of SINR). Let P M {l) be given by 



Pm{i) 



= Ms 7--f e r g (so)+-ferg(0) + 



'l("0) + "2(°0) 



(9) 



In the above equation, I erg (s) is the ergodic mutual informa- 
tion given by 



lerg(s) = trill 



K 



fc=l 



(10) 



A" 



K 



^N k Hz + p k t k {s)) - ^N k r k {s)t k {s) 



fc=i 



fc=i 



rfc(s) 
t k {s) 



Pk 



z+Pkt k (s) 

1 



tr 



K 



(ID 

-1 



Rfc (-lAf +PosR + ^R ? r g (s)J (12j) 



for k = 1, . . . , K. The parameter z takes the value z — 1 for 
the case of the MMSE SINR and the value z — for the case 
of the ZF SINR. The variable so in (0 is evaluated through 
the saddle-point equation 
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1 

M 



I'erq( s o) 



Pn — tr 
M 



(13) 



Ro Im + PoSqRo + ^2 Rfc r fe(«o) 



fe=i 



I' erg (s) is the derivative of I er g(s) with respect to s. The 
expressions of the 0(1) terms vi(sq) and V2(sq) are given in 
Appendix \E\ I erg (0) is obtained by setting s = in I erg (s), 
tk(s), r k (s) above. 

Then for every 7 > 0, the probability density converges 
weakly to Pm{i) in the sense that 



lim M Pm (7) 

M-s-oo 1 



Pm(i)\ < 00 



(14) 



Proof: See Appendix iBl ■ 

Remark 3. As it will become clear in the appendix, this result 
means that for large M the PDF of the SINR becomes asymp- 
totically equal with Pm{i), up to corrections of 0(1/M). 
Remark 4. The solution of ( fTTI ). ( fl2l i has been shown to be 
unique for the case of the MMSE SINR (z = 1) 11311, 11341, 
|35l . To show that this is also the case for the z —> + limit, we 
observe that we can rewrite (fTTTi as r k (s) — (t k (s) + z/pk)^ 1 . 
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PDF of MMSE SINR; N=M=2; P =P,=1; 10 7 runs CDF of MMSE SINR; N=M=2; Pn=P,=1; 10 7 runs 




r(dB) T(dB) 

(a) PDF of MMSE SINR (dB) (b) CDF of MMSE SINR (dB) 

Fig. 1. (a) Probability density (PDF) curves and (b) Outage probability (CDF) curves for the MMSE SINR in dB for N = M = 2, with transmission 
and interference powers po = 2 and p\ = 1 respectively. The channel elements are assumed to be iid. We plot the PDF for the Monte Carlo-generated 
simulations (MC), the LD approximation, the Gaussian approximation and the generalized Gamma approximations. In the Gaussian curves we have used 
NildB.ergi&dB) w ' tn 1dB erg = 10 * log 10 -y e rg and (j 2 dB = o"e r9 /(7er 9 lnlO) 2 . The ergodic mean 7 er9 and variance a^ rg of the SINR can be 
calculated directly, see e.g. [24 1 . The generalized gamma curves have been plotted using the parameters of the generalized gamma distribution as calculated 

in da. 



Hence we may view the ZF z —> limit as the MMSE solution 
with p k -> oo limit (k = 1, . . . K). Since the MMSE (z = 1) 
solution for the rk(s), tfc(s) is continuous with respect the 
values of pk, we may take the MMSE rk(s), tjt(s) solutions 
in the limit of large pk, and then plug them in ( flOl ). setting 
also z = 0. 

Also, note that the most probable value of 7 corresponds 
to the solution of ( fT3l ) for sq = 0. This involves the joint 
solution of ( fl2b . ( fTTT i. which gives the correct value of the 
ergodic SINR El, EH. Expanding the leading O(M) term 
in the exponent of the PDF (i.e. the first three terms) to second 
order in (7 — 7 er g) provides the Gaussian approximation of the 
PDF of the SINR. Furthermore, since Pm(i) in ® is valid 
for all positive 7, not necessarily close to the ergodic value, 
it can provide the tails of the distribution accurately. 

In ||37l we derived a simplified expression for the case when 
all correlation matrices are identical. This result can also be 
obtained from the above analysis by setting N\ = K and all 
other N k = 0, and R = Ri = R, while g ~ CAf(0,l/K), 
Gi ~ CAf(0,diag(pi, . . . ,pk)/K). 

Corollary 1. H37V Let each of the K + 1 transmitters have 
a single antenna with same correlation matrix at the receiver 
given by R. Then in the limit of large K, and N with q = 
K/N fixed the expressions AlOh t \ll\) . t \12\) are simplified to 

I erg (s) = trln[I M +R(>os + r(s))] (15) 

K 

+ ^2ln(z+Pkt(s))-Kr(s)t(s) 

k=l 

K f-J z+p k t(s) 

t{s) = K^l + Rifas + ris)) (17) 

where R it (i — 1, . . . , M) are the eigenvalues of the matrix 
R, while z takes the value z — 1 for the MMSE case and 



z = for the ZF case. As a result, ( 17 3D simply becomes 

7 = Pot(s ) (18) 

with the corresponding expressions for / er9 (0), r(0), t(0) re- 
sulting from setting s — to ( Ii5| ), ( Ii6| ), ( Ii7l ), respectively. The 
expressions of vi(sq), «2(so) ® r e olso accordingly simplified 
( see Appendix \E}. 

To be able to compare the obtained distribution of the 
MMSE SINR with other proposed distributions |[T8), |fl9l , 
[23 1, [24 1, it is instructive to further simplify the assumptions. 
In particular, we have the following 

Corollary 2. In the case of equal power transmit antennas 
Po — Pk — P and uncorrelated receiver antennas R = \m> the 
result simplifies and, to leading order in N takes the following 
simple form: 

P M (7) oc e- K -ilr 7 (19) 

[z+if 

This extremely simple result is quite remarkable. Although 
for large M and close to the ergodic value of 7 this equation 
will behave approximately as a normal distribution, for general 
values of 7 this is far from a Gaussian or generalized Gamma- 
distribution. This is partly the reason why all efforts to 
approximate the distribution of 7 using a central limit theorem 
approach have largely failed, at least for relatively small values 
of M. At the same time, when z = 0, the above distribution 
becomes exactly a Gamma distribution as shown in |20|. 

B. Outage Distribution of 7 

Using the expressions of the probability density Pm(i) 
from the previous section we may now evaluate the asymp- 
totic expression of the outage probability of the SINR 
PM,out{la) = P(7(H) < 70). It turns out that it can be 
evaluated using the information obtained thus far. In particular, 
we have 
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(a) PDF of ZF SINR (dB) (b) CDF of ZF SINR (dB) 



Fig. 2. (a) Probability density (PDF) curves and (b) Outage probability (CDF) curves for the ZF SINR in dB for N = M = 2, with signal and interference 
powers po = 2, pi = 1, respectively. The angles of arrival of the signal and interference paths are 8o = 0° and 8 = 45°, respectively, measured from the 
vertical of the receive antenna array. The angle-spreads of both paths at the receiver are a as = 30°. As in Fig. 1, the LD curve using this approximation is 
consistently closer to the Monte-Carlo generated curve (MC). The way the Gaussian curve of generated is identical to Fig. Hill 



Proposition 2 (Outage Probability). Let Pu,out{l) be given 
by 



(20) 



for 7 < lerg 



and 

1 - Pal out (7) = 



Q (vfeo>o)|*8 



(21) 



when 7 > 7 erg . I' e ' r g( s o) is the second derivative of 
lerg(s) with respect to s. Q{x) is defined as Q(x) = 
J°° dxeT x / 2 /\/2n. The definitions of sq, I e rg(so) and V\{sq) 
can be found in Proposition Q] and Appendix \E\ The de- 
pendence of sq on 7 can be obtained through U3h ^erg 
corresponds to the value of 7 in ( 1731 ) when s — 0. The 
parameter z = 1 for the MMSE (0 case and z — + for 
the ZF (0) case. 

Then for every 7 > 0, the outage probability function 
converges to PM,out{l) ' n the sense that 



M 



lim M |P M (7(H) < 7) - PM,out(j)\ < 00 



(22) 



Proof: See Appendix [C] 



modulations can be expressed as 

Q{V2j) BPSK 

p e ( 7 ) = { QiV^l Q p SK 



log 2 L 



QW-Bi^ L-QAM 



(23) 



where the latter expression holds approximately for large L 
J25). The average BER is given by the following 

Proposition 3 (Average BER). Define the following function 

B ERm = ^e I "' 9 ^~ I ' r9 (^ )+ ^ vl (^ )+ ^ jI '< !r 3 < - 217(^24) 

. r ( —r (—• 

where T(x,y) = dtt x ~ 1 e~ t /T(x) is the normalized in- 
complete Y -function and I erg (s), P erg (s) and «i(s) are the 
ergodic mutual information, its derivative with respect to s, 
and the variance defined in MOh ( liiD and ( 1401 ) respectively. 
The above function and parameters are defined both for the 
MMSE (z = 1) and the ZF z — 0) receiver cases. Also the 
parameters a, b describing the modulation are defined in ( 1251 ). 

Then if BERm is the average uncoded bit-error rate of the 
MMSE (0 and the ZF (0) receivers, in the limit of large M 
we have 



lim M BERm 

M— >oo 



BER M \ < 00 



(25) 



Proof: See Appendix ID] 



C. Evaluation of Average BER 

In addition to the outage probability, another important 
metric of performance for the linear receivers is the average 
uncoded bit-error probability (BER). This can be expressed 
as an average over 7 of P e (7), the bit-error probability 
conditioned on the channel realization, which for different 



IV. Numerical Simulations 

To test the applicability of this approach, we have performed 
a series of numerical simulations and have compared our large 
deviations (LD) approach with Monte Carlo (MC) simula- 
tions, the Gaussian approximation and the generalized gamma 
approximation by [23|-|26|. We start with the simpler case 
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where no correlations are present in the receiver side using 
different powers for the transmit antennas. In Fig. 1 we plot the 
probability density (PDF) and the outage probability (CDF) of 
the MMSE SINR in dB for the 2 x 2 antenna case. The PDF 
curve of our large deviations (LD) approach is consistently 
closer to the Monte-Carlo (MC) numerical curves. The same 
is true also for the outage curves even for such small antenna 
arrays. 

In Fig. 2 we plot the PDF and CDF curves for the zero- 
forcing (ZF) SINR in dB for the 2x2 antenna case, using 
different correlation matrices R for the two transmitter paths. 
In particular, we parameterize the correlation matrix elements 
using the mean angle of arrival 8, as measured from the 
vertical of the antenna array, and a Gaussian angle-spread a as 
as follows: 



-Rah — C 



d(j)e 



2-iridab sin(0) /A 



(26) 



where A is the carrier wavelength, d a b is the distance between 
antennas a,b, taken to be d a f, = (a — b)X/2 and C a 
normalization to ensure R aa = 1. Using the above notation, 
the angles of arrival of the signal and the interferer are 8q = 0° 
and 9 = 45°, respectively, while all angle spreads are taken to 
be Gas = 30°. In this case, we also see very good agreement 
with the Monte-Carlo curves. 

Finally, in Fig. 3 we test our predictions of the uncoded 
BER, both for MMSE and ZF. In Fig. 3(a) we take un- 
corrected receivers and compare to Monte-Carlo simulations 
and the generalized gamma approximation. We see that at 
large SNRs, the generalized gamma distribution deviates up 
to several dB. In contrast, our LD approximation is quite 
close to the numerical curve. We see similar behavior for our 
approximation in the ZF case. In Fig 3(b) we plot the BER 
as a function of angle-of arrival of the signal path, in the 
presence of two interfering paths, for several angle-spreads 
and receive array sizes. We find that low angle-spreads lead 
to deterioration of the BER when the signal path has the same 
direction of arrival as the interfering paths. In addition, we 
find that lower angle-spreads increase the BER away from the 
interferers' direction. This last observation is due to the fact 
that higher angle-spread leads to higher diversity and hence re- 
duced outage probability. Interestingly, an angle-spread of just 
cr as = 5° is enough to make two interference paths separated 
by 30° practically indistinguishable for M = 6, N = 3. 

V. Conclusion 

In this paper we have used a large deviation approach to 
calculate the key statistics of the SINR, i.e. PDF, outage 
probability and BER for the MMSE and ZF receivers of 
the Gaussian MIMO channel with arbitrary receive antenna 
correlations. Our results agree very well with simulations both 
close to the peak of the distribution as well as at its tails, where 
other suggested approximations, such as the Gaussian or the 
generalized Gamma distributions are inaccurate. As a technical 
byproduct, we have found an exact relationship between the 
SINR distribution and the moment generating function of a 
difference of related mutual informations. Remarkably, the 
accuracy of the calculated distribution, even at its tails, is a 



by-product of the robustness of the Gaussian behavior of the 
MIMO mutual information. Several direct generalizations are 
possible. This approach may be generalized to include multi- 
tap or frequency selective MIMO channels 

Appendix A 
Proof of LemmaQ] 



The moment generating function of 7 is 



9m{s) = Eh 
We can integrate over go to obtain 



(27) 



9m(s) 



E 



H11 



dg e 



-gj(lM+spoRS /2 L _1 Rj /2 )go 



-In dot ri JU +p sRo /2 L- 1 Ry 2 l 



(28) 



<: Ho {e- A ™} 



where the quantity AI(s, Ho) is exactly (jTj. AI(s, Ho) 
and therefore gu (s) will be analytic in s when 3t(s) > 
— A m ;„(Ho), where A m i„(Ho) is the minimum eigenvalue of 



-1/2 



We will assume 



the matrix p^Rg 1/2 (l M + H H J R ( 
that in the large M limit A m i„(Ho) will converge with proba 
bility one to a fixed value \* nin and hence for 5ft (s) > — A^ in 
gM^oa(s) is analytic. This has been shown for K = 1 
and is expected to be true for general K > 1 



Appendix B 
Proof of PropositionQ] 

Before discussing some elements of the proof, we introduce 
the normalized mutual information difference as 

1 



«( S ,H ) = — A/(s,H ) 



(29) 



where AI(s, Ho) is given by (0. We also introduce an 
important property of SI(s, Ho). 

Lemma 2 (Hardening of SI(s, Ho)). In the limit M — > 00 
the quantity <5/(s, Hq) converges with high probability to 



5I erg (s) 



J-erg (^) ^erg (0) 



M 



M 



(30) 



where I erg (s) is defined in ( |70| ), ( 1771 ), ( 1721 ). 



This Lemma was proved in ||34l for the case s £ M + . We 
will assume it is valid for 3t(s) > — A* liTl . We should mention 
that for the case K = 1, or for the case of equal correlation 
matrices Rfe = Ro (k = 1,...,K), the generalization to 
5ft(s) < can be inferred from [39 |. From the above result and 
using the linearity of the derivative operation, we can deduce 
the "hardening"of all derivatives of SI(s, Hq) with respect to 



Corollary 3. (Hardening of Derivatives of 81{s, Ho)) In the 
limit M — > 00 the derivatives of 51{s, H ) with respect to s 
converge with high probability to their deterministic equiv- 
alents, which are the corresponding derivatives of I erg (s), 
defined in ([70}, (ED, GUI. 



7 



From the convexity of the function AJ(s, Ho) with respect 
to s, we can deduce that I'J rg (s) < 0. 

To show Proposition Q] we start by expressing the probabil- 
ity density function of 7 as follows 



P M {l) 



M 



ds 
2it~i 



; Ms 7 



9m{s) 



by fk) are functions of the realization of H , directly or 
through s i which is the solution of dl3t . Nevertheless, from 
Corollary [3] we can replace the derivatives of 5I(so,Ho) with 
their deterministic equivalents to leading order. As a result, to 
leading order in M we have 



(31) 



Pm{i) 



Me So7 



E 



Hn 



M 



ds 



where 



f(s,K Q ) = sj-6I(s,H ) 



(32) 
(33) 

Keeping in mind that in the large M limit 5I(s, Ho) =0(1) 
we proceed to first integrate over s before averaging over Ho. 
Since for 5ft(s) > — A m j n (Ho) f(s) is analytic, we deform the 
contour of the s-integral to pass through the saddle point(s) 
of f(s) from the steepest descent path [40|, which are defined 
by f(s ) = or 



^2* |E [<$/"(*„) 



H„ 



-M57(so,H ) 



(38) 



y/2TT\E[SI"( S0 )]\ 



g(so) 



7 = 5 1' (so) 



Po , 
—— tr 

M 



Rn I 



LAf 



P0S0R0 + H H ;z 



(34) 



It is easy to see that the above equation only has real solutions 
for Re(so) > — A m j„(Ho). This is so, because in this region 
the right-hand-side is real only if sq is real. Also, since the 
right-hand-side above is a decreasing function of s, (becoming 
unbounded when s — > — A m j n (Ho) and going to zero when 
s — > 00), it can also be shown that it can only have one 
solution, which depends on 7. Hence for M —> 00 we expect 
the resulting limiting equation 7 = 5I' erg (so) to have a single 
real solution for s > — A* rain . 

For large M, the integral will dominated by the behavior 
close to the saddle point. As a result, we may expand the 
exponent close to sq. Thus 



To conclude the calculation, we need an expression of g(so). 
Clearly, the "hardening" of the mutual information itself 
SI(so, Ho) has also been shown elsewhere. However, here we 
need an expression accurate to 0(1/ M), hence we will need 
the next, i.e. 0(1) correction. This correction can be evaluated 
using the fact that A/(so,Ho) is a difference of two MIMO 
mutual information functions with noise covariance matrix that 
differs by soPoRo- We can then take advantage of a number 
of works in the literature that has analyzed the statistics of 
mutual information functions. 

Lemma 3 (CLT for AI(s, H )). In the limit M -> +00, 
N k -> +00 (for k = 1,...,K), such that n k = N k /M 
remains finite, and for s > — A* lin the quantity A7(s, Ho) 
in becomes asymptotically normal. In particular, 

A7(s,H ) — AI erg (s) 



(39) 



m = /( So )+E(- 1 ) fe ( s - s o) fe /fc 

k=2 



(35) 



fk 



M 



PqSqRo 



HoH^- 1 



VM s ) 

where AI erg (so) and its related parameters are given by ( liOl ). 
The variance V\(s) of AI erg (s) is given by 

wi(s) = var(7(s,Ho))+vax(/(0,H )) (40) 

- 2cov(/(s,H fe ),/(0,Ho)) 
= -logdetllx-naSal 

- logdet|I K -n So| 
+ 21ogdet|Ix - niSi| 

The elements of the positive-definite matrices n and £ are 
given below 

2 



Since fi < the steepest descent path in the neighborhood 
of so is s = sq + it, tel. Keeping the first non-trivial term 
in the expansion of f(s) d36l ) in the exponent, we expand the 
rest obtaining an expansion of the form 



M 



Mf(so) 



dte 



1 + y Mi A q (t)\ (36) 



9=1 



where the function A q (t) can be expressed as an expansion 
of t, with the minimum degree 3q if q is even and minimum 
degree 3q + 1 if q is odd. Integrating over t and performing 
simple power counting of M we conclude that to leading order 
in M we have 





= Sab 
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So,ab 


= Sab 
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Ho, ab 






= S a b- 
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ni ia b 


Aa 



Pa 



Z +Pata(so) 



Pa 



(41) 



Z + Pat 0,(0), 

tr [R^Q^-^QCO)- 1 ] 



Pa 



(Z + p a t a (so))(l + Pata(0)) 



Pin) 




3 M( S0 7-<57( So ,H )) 



1 



0(^; 



for a, b = 1, . . . , K and z = 1 (z = 0) for the MMSE (ZF) 
cases. The matrix Q(s) is defined as 



A number of comments are due for this expression. First, at 
least in principle, 6I(sq,Ho) and all its derivatives (given 



Q(s) = I 



M 



-Po 



sR + y^Rfc?-fc(s) 



(42) 



For convenience we generalize the above notation to include 
^-2.ab, when any of its indices a, b can take the value 0, in 
which case the corresponding matrix R a ( and/or R;, ) becomes 
Ro and N a= o — ¥ M. 

Although we do not formally prove this lemma, we will 
briefly motivate its validity and discuss how one can go about 
to prove it. The Gaussian behavior of MIMO mutual informa- 
tion functions was first introduced in [31], where in addition 
to the ergodic mutual information of the form appearing here, 
the variance of the difference of two mutual informations in 
both of which the same random matrix appears was calculated 
using the replica trick with both complex and Grasmann vari- 
ables. Using this methodology the variance t>i(so) above was 
evaluated. Furthermore it was shown that all higher cumulant 
moments vanish as increasing inverse powers of M _1 . This 
shows that AI(s, Ho) converges to a Gaussian variable in the 
large M limit. Similar results have been shown using more 
formal arguments for the case of a single mutual information 
function with Kronecker-correlated Gaussian channels in [ 32 1 
or with independent but not identically distributed channels 



Armed with the above result, we can now integrate over the 
channel Ho by changing variables, from Ho to the random 
Gaussian variable Z — (A/(so,H ) — AI erg (so))/ ^ vi(s ). 
The reason we shift from A/(so,Ho) to Z is because we 
know from the analysis above that it is Z that becomes 
asymptotically Gaussian. It is also the case that gAi(s) itself 
involves the expectation of an exponentially small quantity 
(e~ A/< - So ' H °)) when M is large, hence its average is not 
necessarily well definecH 



g(s )e 



AI crg (s ) 



E H , 



-(AZ( So ,H )-A/er 9 (so)) 



(43) 



3 -Z«i(s ) 



dz 



■V- *«i(»o) (1 + 0(M" 1 )) 



e ^(l + 0(M- 1 )) 



The corrections of order 0(1/M) stem from a number of 
sources. Specifically, the correction to AJ(so, Ho)— AI erg (so) 
is 0(1/ M) ||3~n , while the correction to the variance i>i(so) is 
0(1/ M 2 ) EQ. Both these corrections result to an 0(1/ M) 
correction to the above result. Also, for finite large M we 
may incorporate corrections to the Gaussian approximation by 
including the higher order statistics, e.g. the skewness [41 1. 
Here again the leading contribution stems from the skewness, 
which is 0(1/ M) IT3T1 . 

The expressions in Lemma [3] allow us to express the second 



2 For example, even if Z above is asymptotically Gaussian the expectation 
E [e~ MZ ] is not well denned. 



derivative of 8I erg (so) with respect to s as follows: 



Hernia) 



-po n 2,o, - Pon 2 ,oo (44) 

3=1 
K 



2,00 



-v 2 (s ) 



The second equality follows from the expression of the deriva- 
tives of rk(s), tk(s) in dTTb . ( TT2b with respect to s in terms of 
the matrices U2, £2- Finally, the last line above defines V2(s) 
in ©. 

Appendix C 
Proof of Proposition!!] 

We will now provide some details in the proof of <|20t . We 
will deal only with the case 7 < -f erg , since the opposite 
case 7 > -f er g can be analyzed in a similar way. PM,out("l) is 
defined as 



PM,out(l) 



dtP M (t) 



(45) 



up to negligible corrections 0(M^ 1 ) due to replacing Pm(4) 
for Pjvf(i). The analysis is based on the fact that for large M 
the outage probability PM,out(l) is determined from the be- 
havior of Pm (t) close to the We will need to focus separately 
in two regions of interest in the interval 7 € (0,7 ers ]. In the 
first region I7 — 7 erg | — 0(1), to asymptotically evaluate the 
outage probability we expand the exponent of Pm(^) ® m 
7 around the end point of the integral. Since Pm(i) is an 
increasing function for 7 < j erg its derivative will be always 
positive in this region. Hence we have 



P 



M.out 



(7) 



dt 



o Ms (t-j) 



V2tt 

, e M SO 7-AJ CT! ,( S o)+ ' ll '" ) l" 2( ° 01 
M So j-AI srg (s ) + ^^p2(fol 



(46) 



2ttMs 



(l + 0(M- 1 )) 



where so above is evaluated at the endpoint 7. We have 
used the fact that to leading the derivative of the exponent 
with respect to 7 is simply Mso. The above approximation 
begins to break down when |s |\/M <C 1, i.e. in the region 
>/M|7 — 7er ff | = 0(1). Although this situation will rarely 
occur when we take the limit M — > 00 for fixed < 7 < 7 erg 
it useful to pay attention to this region so that we can provide 
an approximation that is valid for every 7 when M is large 
but fixed. In this increasingly diminishing region as M — > 00, 
the Gaussian approximation of the SINR is valid, where the 
probability density of Pm(^i) will be approximately quadratic 
in 7. Hence we expand the exponent of Pm (t) to second order 
around the endpoint t = 7, and then integrate over t. After 
some algebra and using the fact that 



ds (j) 
c?7 



1 



<tf£>o(7)) 



(47) 
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MMSE BER; N=M=3; 10° runs 



MC 

LD 

- - Hong et a 



ZF BER; N=3:0 , -[30 GC\ I 




SNR (dB) 



(a) BER for MMSE SINR (b) BER for ZF SINR 

Fig. 3. Uncoded BER for (a) MMSE and (b) ZF SINR. In the MMSE case, we take M = N = 3 and compare to Monte Carlo simulations and generalized 
gamma distributions with iid channels. In the ZF case, we take N = 3 and M = 6 or M = 9. We plot the BER as a function of the angle of the signal path 
in the presence of two interfering paths arriving at 8 = 30, 60°, with various angle spreads cr a s = 1°, 3°, 5°. The dashed lines in (b) are the Monte Carlo 
simulation results. 



we obtain d20} , To obtain the expression in d2lT i we express 
Pout (70) = 1 — P(7 > 70) and work as above with P(7 > 70)- 
The final expressions ( f2Qb . d2lT i smoothly interpolate between 
(ffoT i (for I7 — jerg\ = 0(1)) and the Gaussian approximation 



(for 1 7 - 7 e 



0(1)). 



Appendix D 
Proof of Proposition^ 

The average uncoded bit-error rate (BER) for signals with 
modulation as in d23l l can be expressed in terms of the 
moment-generating function g(is) as follows 



BER 



M 



b I d 7 P M (7)Q(^7) 
/o 

b 1 r+00 




(48) 



ds g M (is) 



, 2ni s + i0+ - i2sM/a J 

dt 9M ( 2M~ ) 

, ~7t t + i 

In the first equation the parameters a, b correspond to the 
different cases in (f23t . The second equation results from 
the definition of Pm(i) in terms of the moment-generating 
function. The third equation follows by deforming the integral 
from the real axis to follow the branch cut appearing due to the 
square root. Using d43b to express <7m(s) in terms of AI erg (s) 
etc, we get 

b f°° dt 

BER 



- e - A/ "s(w)+2 t, i(w ) 



dt 



(49) 



2tt J y/t{l + t) 
b 

•OO 7 j 

. e -2irC s (OT)*(i + 0(M -1 )) 



lo Vt{l + t) 

In the second line above we have expanded the exponent for 
small arguments and kept only the 0(M) and 0(1), neglecting 



all lower order terms. Integrating the above expression over t 
gives (124-b . It should be noted that if we wanted to be strict 
regarding the leading corrections being 0(M _1 ), in the above 
expression the arguments of v\{s), I erg (s) and I' erg {s) should 
be set to s — 0, rather than s = a/(2M). Nevertheless, 
we have found numerically that these expressions are slightly 
more accurate. 
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